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We examine dynamical mass generation in QCD with large current mass quarks. 
A renormalization group analysis is performed to separate fermion self-mass into a 
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dynamical and a kinematical part. It is shown that the energy scale of the Schwinger- 

i-G ■ 

Dyson (SD) equation and the effective gauge coupling are fixed by the current mass. 
• i— i . 

The dynamical self-mass satisfies a homogeneous SD equation which has a trivial so- 
lution when the current mass exceeds a critical value. We therefore suggest that the 
quark condensate, as the function of the current mass, observes a local minimum around 
eA QC D- 
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1 Introduction 

Dynamical chiral symmetry breaking in QCD has been extensively studied [|TJ. The standard 
tool to study this problem is the Schwinger- Dyson (SD) equation, i. e. the fermion gap 
equation. Since one is usually interested in the chiral limit one studies the gap equation in 
the limit of the vanishing current quark mass. In the presence of a quark whose current mass 
is much larger than Aqcd, the nature of solutions to SD equation may exhibit non- analytic 
behavior in m. It has long been argued by Pagels || that if the current mass is large, it is 
possible that m ^ does not belong to an analytic extension of m — 0. In this article we 
explore this scenario using the renormalization group equation (RGE) analysis. 

Recently, Langfeld, Alkofer and Reinhardt || have reported on a numerical study of the 
Schwinger-Dyson equation for massive quarks in the background field of a classical vacuum 
solution. Besides the usual spontaneous chiral symmetry breaking solution for m = 0, they 
find that the total quark condensate as a function of the current mass exhibits a discontinuous 
drop at m c ~ 70MeV. This is interesting because not much has been known about the 
behavior of quark condensates for quarks of large current mass. 

However, in a recent conference report Q we have discussed the problem of dynamical 
mass generation in QCD with large current mass (m ^> Aq C d)- It was pointed out that using 
the renormalization group analysis one can separate unambiguously the quark self-mass (the 
part of the self-energy which commutes with 75) into a dynamical and a kinematical part. 
The dynamical part Bp satisfies a homogeneous SD equation but with an effective coupling 
constant g 2 (t) ~ (In A m ) -1 and hence when m ^> Aq C d it describes weak coupling regime. 
As a consequence, the SD equation for B^ can be solved quite reliably in a single effective 
gluon exchange and one finds that Bp has only a trivial solution when m exceeds a critical 
value m c which is of the order of 2.7Aqqd- The kinematical part Bx, however, is a non- 
singular function of the current mass and satisfies an inhomogeneous equation and hence 
is proportional to m. The total fermion condensate can be written similarly as a sum of a 



dynamical and kinematical part which is expected to take a drop near current mass of the 
order of m^ . 

The aim of this paper is to elaborate on our method which utilizes the renormalization 
group techniques. The dynamical self-mass may be understood as arising from an effective 
interaction with the gluons. It would be of interest to explore if the dip in condensate has 
any observable consequences ||. 

2 The Renormalized SD Equation 

The quark self-mass B{p 2 ) satisfies the following SD equation (the gap equation) in the 
momentum space 

B(P 2 ) =m+ ig2C f N) Tr J tS^Cp " k) 7f J5(k)T v {p, k) (1) 

where S(k), D fJ,v (k) and T u (p, k) are the complete quark, gluon propagators and the proper 
vertex respectively, m is the current mass; and C^{N) = J ^~ 1 for SU(N). The quark 
propagator S is related to B by 

S-\k)=t/+#{k)-B{k 2 ) (2) 

where $(k) is the wave function correction factor. Before trying to solve ([!]), we must rewrite 
it in terms of renormalized quantities and specify a renormalization prescription. We take 
the point of view that the renormalization constants in the renormalized SD equation may be 
defined as in perturbation theory. We also require that after carrying out the renormalization 
prescription we get a finite SD equation. We adopt the dimensional regularization and mass 
independent renormalization scheme 11. We define renormalized quantities as follows 

S(p;g,m,£;e) = Z F (g(fj,)\e)S R (p\g(n),m(iJ,),€(iJ,);ti); 



1 We have neglected the non-trivial topological gauge configurations. Inclusion of them may lead to a 
contribution from the gluon condensate. However, the essential feature of the behavior of the total quark 
condensate should not be substantially affected by this. 



D^ u (p;g } m } £;e) = Z A (g(fi); e)D%(p; g(ji),m(ji),£(ji); //); 
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F»(p;g,m,£;e) = Z F Z A T li {p;g(fi),m(fj,),^(fiy,fi) (3) 

B(p;g,m } £;e) = Zp X B R {p, g(jt),m(ji),€(ji); //) 

where 

#(e) = Z g (g(ii);e)g(n); 

m(e) = Z m (#(/i);e)m(/i); 

£(e) = Z^O*);e)^). 

The renormalized quantities are functions of g(/i) and e. £ in (§) is the gauge parameter. 
Substituting @ into (HD we obtain the renormalized integral equation (in Euclidean 4 — e 
dimensions) 

E'W) = Z F (e)[Z m (e)m^) + Z 2 g (e)Z 1 J 2 (e)g 2 ^)C 2 (N) (4) 



:Tr / {2ttY-^ D ^ <J) ~ k] ^ Sr ^ ; ^ T R^ k] ^)1- 



4-6 7 (2tt)< 

In ([|) we have written S#(p; g(/i), m(/x),£(/z); /x) as Sr(p; /i) for short. 

Eq. dj) is in general not self-consistent since its right hand side contains divergent renor- 
malization constants as e — > 0. It has been suggested by Johnson, Baker and Willey that 
a suitable choice of the gauge parameter £ (Landau gauge) would lead to a finite Zp as e — > 0. 
Subsequently, it would be possible to require that the divergences arising from the integral 
cancel those in Z m (/j,, e)m(/i) and make (§) finite. Thus, a consistent solution to (§) should 
satisfy the condition of finiteness of Zp(fi, e) as e — > 0. We show below using RGE analysis 
that when the quark current mass is large, this condition can be approximately satisfied. 



3 Renormalization Group Analysis 

Equation (0) contains the exact propagators and vertex. If the current mass is small, 
B R (p; g, m, £, \i, ) may be expanded in powers of m(/x) 

B R (p;g,m,^fx) = B R (p; g, 0, £; //) (5) 

+ m(ii)B?(jp; g, 0, £; /i) + m(/x) 2 5^(p; #, 0, £;//) + •• • 

where we call the first term B Q as the 'dynamical' part of B R and other terms as the 
'kinematical' part. A non-vanishing dynamical part signals the spontaneous chiral symmetry 
breaking. @ is just the chiral perturbation expansion, which can only make sense when the 
series is convergent or m is small (m <C Aqcd)- We would like to develop a scheme applicable 
when m > Aqcd, when the power series solution does not converge. A convenient method 
in this case is to apply RGE treatment. We may write (Q) as a functional equation 

F(p; g(fi), m(fi), £(//); /i) = 0. (6) 

From (|J) it follows that 

(l*S- + MS- + 7mm/- + 5(g)£)F = 0. (7) 

on og am o^ 

In a mass-independent renormalization scheme, /3(g), j(g), and 5(g) depend only on g(p). 
Introducing t by 

i = ln™* (8) 

n 

with map = m(\x = rap) and observing that F is homogeneous of order 1 in p, m and //, we 
have 

{ ¥t +p T v + ^~ l)F(p; 9 > eW ' e; /i) = °- (9) 



In (|^) m'(/x) = e t m(fi). From (0) and (|9|), it follows that F satisfies the following RGE: 
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+ 6(g)^ + l)F(p;g,e t m',^ f i) = 0. (10) 

From (|Ty) it follows that 

F(p; (y, eW, £; /i) = e^e"'?; g(t)M(t),Z(t);n) (11) 

where the effective parameters g(t), m'(t), and £(£) are defined by 

m , 

f dx ,.,„ s 

t= jm (12) 

m' (t) — mf Y/i) exp / dx ™ , (13) 

7 p(x) 



and 



It should be noted that even though £ is a function of /x, the effective coupling g(t) depends 



on mp only. An explicit calculation of (|12|)-(|I1) gives the following results in one-loop 



approximation in the minimum subtraction scheme 

9 2 it) 2n 



a ( m P) = ~^r = a i„ m P , (15) 
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m 



and 



m'(t) = e m(u) exp / dx = e m(t) = e mi^i = m P ) = /i, (16) 

J p(x) 



*(*) = 1 " TjrA^—r, (17) 



Colin 
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where (3q = ~ n/ , d^ = 2(33 _ 2 " / , [|J and £o is a constant of integration. It is interesting to 
observe that m'{t) is just the scale parameter /i. This arises because we have defined t by 
the on- shell current mass irtp. 

The physical consequence of the RGE analysis follows from equation (|ll]). Instead of 
solving the SD equation ([J) with a large current mass m directly by approximation procedures 
for r„ and D^ v , we can solve, equivalently an effective SD equation govern by a much smaller 
mass fi (= m'(t)) and running coupling and gauge parameters. A suitable approximation to 
these are given in fll5|)-([T^). It is seen from fllTD that when mp ^> Aqcd, £{t) — 1 (Landau 



gauge) and g(t) — > as m P — ► oo. Thus a large current mass justifies the single gluon 
exchange approximation for the effective SD equation and ensures the compatibility of the 
condition Zp{g, f, e) — ► 1 + 0(g 2 ) as £ — > 1, whereas, it would have been hard to justify such 
an approximation in ([!]). 

As a consequence of these results, the renormalized Green's functions D^ and T^ R may 
be approximated by 

Dr(p; 9, »,€;») = 5 + °Gr) 

r^(p^,/i,e;/i) = 7 M + o(^ 2 ). (is) 

Substituting (|T8| ) into the effective SD equation as given in (|Tl|), we arrive at the RGE 
improved effective SD equation to 0(g 2 ) 

B (fl = *.(* <)„ + itfft W / g i^l^JgL + o^), (19) 

where we have used the notation B(p 2 ) = B R (p;g,fj,,^;fi). Finally, from a solution to (|Kp, 
we can reconstruct B R (p; g(/x), m(/z), £(/z); /J,) defined in Eq. (|j) from the following equation, 
which is easily derived 

B R (p; g(ji), e 4 m(/i),e(/i);/i) (20) 

= e'exp[- J dx^}B R (e- t p;g(t),^m^)- 



Now in ( |T9| ) the inhomogeneous term contains a mass of the order of Aqcd, we can thus 
use a chiral perturbation theory to expand B R (p;g(t),mo,^(t)]fj,) as power series in m Q (in 
our case, of course, m^ = /i. We use a different symbol to make this decomposition more 
obvious) 

B R (p;g(t),m ,Z(t);n) = Bf m ip;g(t),0^(t); p) 

+moB* n {p;g{t),m ,Z{ty,v) (21) 

where we have called the term independent of m as the dynamical self- mass and have lumped 
the rest as the kinematical part. In this way of splitting, it is obvious that the kinematical 
part is a regular function of the m$ and it goes to zero as m$ — > 0. But what is crucial is 
that both Bdyn and B^ in depend on nip in a non-analytical way through g(t). 

Substituting ( pT|) into ( p0|) we obtain the following structure for the self-mass for a heavy 
quark (taking 7^ = 0) 

B R (p;g(^,m(fi) J ^);fi) = e*S d j /n (e~V;y(*)) 

+m(/i)exp(- / dx ™ )B kin (e~ t p;g(t) } m ). (22) 
J g (3{x) 

It is to be emphasized again that the separation of the self-mass into a dynamical and a 
kinematical part is not the usual power series expansion in current mass. In fact, substituting 
( |2l~D into fll9|) we derive an effective homogeneous SD equation for B dyn 

B dyn ( P ; g{t)) = 3f(t)C 2 (N) f £* ^ ,^twL + 0{f) (23) 



(2vr) 4 (p - fc)2 fc2 + £2^2 

while £?fci n satisfies 

**.<*««» - 1 + s-ftw / gi (p r t) , (p fgg^ (*■ - 3^)- (24) 

Both (E3p and (123) have been studied in the context of chiral symmetry breaking [RJ . 



4 Solutions and Discussions 

The effective SD equations for the heavy quark differ from the those for the light quark where 
the chiral symmetry breaking is of concern. In the latter case the gauge coupling constant is 
a function of the momentum transfer. In the infrared range the coupling becomes arbitrarily 
large and a non-trivial solution exists for light quarks. In (|23|) and fl24|) the coupling constant 
is fixed and allows us to seek for a solution in the wider range (not just the asymptotic solution 
in the limit p 2 — > oo). However, a non-trival solution to (|23|) should not be referred to as 
the signal of chiral phase transition. 

Eq. (^||) is equivalent to the following differential equation 

X A 2 + 2X A + A |R2 M =Q 25 

ax 1 ax x + Bj yn (x) 

together with the boundary conditions (x = p 2 ) 

= (26) 



% 3 r Ody n \Xj 



2 CUD(ly n [X 1 

0; x 



x=A " x=5 



where A —> oo and 5^0 and A = 3g C2(N)/16ir . We are interested in finding not so 
much in existence of non-trivial solution to (|25| ) but in determining when the trivial solution 
Bd yn {x) = is in fact a stable solution. To examine the stability of the trivial solution, let 
us linearize ( p5|) about B dyn = e(x). We find 

x 2 e(x) + 2xe(x) + Xxe(x) = (27) 

with 

xe(x) + e(x) \ X=A = 0; x 2 e(x) [^ = (28) 

whose general solution is 

e{x) = c x x p + + c 2 x p - (29) 

where 

1 1 



p± = --± -Vl^AX. (30) 



c\ and C2 are to be determined by the boundary conditions. If c\ and c% are non-zero, then 
the solution e(x) — > oo as x — > and thus the trivial solution is not stable. Thus stability 
requires c\ = ci = 0. Substituting ( p9|) into the b.c. (|28|), we find the following condition for 
the trivial solution to be stable: 



1 + V1-4A i, v^ , 01 s 
= i — i 2 . (31) 



1- y/l-AX A' 



If A < 1/4, a solution to (31) is not possible since the left hand side of ([H]) is finite while the 



right hand side tends to zero in the limit 5 — > and A — > oo. If A > 0, yl — 4A is purely 
imaginary, then (|31~1) is a transcendental equation for A [[l(J. For fixed A and 5, it has an 



infinite set of solutions for A which becomes dense over the whole domain A > 1/4 as A/8 
becomes large. Hence the critical point is A c = 1/4 or a c = tt/4. For A < A c we have a stable 
trivial solution B^ yn (x) = 0. The critical current mass is 

m c = Aqcd exp - — ~ ek QCD . (32) 

(3 a c 

Implicit in the analysis above is the assumption that there are no other stable solutions. 

The solution to (|24|) for B^ in can be derived by the iteration process while Bd yn = must 
be substituted into fl24|) . To the first order one has Bki n = 1 + 0(g 2 (t)) and the kinematical 



part of total self-mass reads from (22 ) 



B* (p;g(ji),m(ji),Z(n);n) = m(//)(6 In -)- c / ft (l + 0(f)) (m P » A QCD ) (33) 

/i 

where c = 3C 2 (iV)/87r 2 and b = (3 /8ti 2 . 

We have concerned ourselves with the self-mass for the heavy quark. Having assumed the 
chiral phase transition for light quarks, we may draw some conclusions for the self-mass of the 
massive quark as the function of its current mass. In the range of m ~ 0, the B^ dominates. 
As m becomes large, the kinematical part grows almost linearly and the dynamical part 
changes slowly. At mp = m c , the dynamical part drops to zero and the total self-mass will 
exhibit a local minimum. As m gets even larger, the self-mass is completely governed by 

10 



the kinematical part (except for the possible contributions from the gluon condensate) and 
eventually blows up as shown in (|33| ) when m — > oo. The quark condensate defined as 

d A k 



{m) = -JT^yTrSFik) 



d*k Tr B« _ r _d^ Tr _ B% 



(2tt) 4 k 2 -(B% + B%) 2 J (2tt) 4 k 2 -(B% + B%) 2 

(qq) D + (qq) k ( 34 ) 



is expected to observe the same drop when Bp = 0. 
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